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Previous authors have calculated the mass function of primordial black holes only on scales which 
are well outside the horizon at the end of inflation. Here we extend the calculation to sub-horizon 
scales, on which the density perturbation never becomes classical. Regarding the formation of black 
holes as a 'measurement' of the (high peaks) of the density perturbation, we estimate a mass function 
by assuming that black holes form as soon as inflation ends, in those rare regions where the Bardeen 
potential exceeds a threshold value of ^ c — 0.5. 
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I. INTRODUCTION 

During inflation the vacuum fluctuation of each light field is promoted around the time of horizon exit to a classical 
perturbation. According to present thinking, one or more of these light field perturbations is responsible for all 
primordial perturbations in the Universe Q . It is usually supposed that primordial perturbations are absent on those 
very small scales which are still inside the horizon at the end of inflation, on the ground that the vacuum fluctuation 
on such scales never becomes classical. 

These ideas are usually applied to the perturbation on cosmological scales, which at least in the standard cosmology 
are far outside the horizon by the end of inflation. A different application though is to the formation of Primordial 
Black Holes (PBH) 0,0. After inflation, as each scale comes back inside the horizon, an appreciable number of 
black holes (BH) will form provided that the mean-square energy density contrast smoothed on that scale is not too 
far below unity. The black holes form in those relatively rare regions where the density contrast is actually of order 
unity, and their mass is of order the total energy inside the horizon. Assuming a specific form for the spectrum of 
the density contrast, this allows one to estimate the mass function of the black holes using for instance the Press- 
Schechter formalism 4] 1 . PBH abundance calculations provide a powerful tool to study inflationary models through 
their primordial power spectrum over a wide range of scales which are difficult to constrain by other data from Large 
Scale Structure surveys or Cosmic Microwave Background experiments 0, ITU . ITU Ip3 | . For a comprehensive recent 
review on PBHs see e.g. |T^ . 

The standard PBH calculation makes sense on those scales which are well outside the horizon by the end of inflation, 
but what is supposed to happen on smaller scales? As the comoving wavenumber k is increased above that wavenumbcr 
fc cn( j which leaves the horizon at the end of inflation, the standard calculation ceases to make sense and a different 
calculation must take its place, which takes on board the essentially quantum nature of the vacuum fluctuation. In 
this paper we explain for the first time how such a calculation should proceed. 

The layout of the paper is as follows: in the next section we briefly review how fluctuations in the inflaton field are 
generated during inflation and calculate the spectrum of the fluctuations at the end of inflation. In Section ITTll we use 
the primordial power spectrum as input for the Bardeen potential and calculate how the latter then evolves after the 
end of inflation during radiation domination. In Section llVI we investigate how black holes are formed on sub- horizon 
scales using the Press-Schechter formalism and calculate the mass function of these black holes. We discuss our results 
and conclude in Section Ivl 

II. THE VACUUM FLUCTUATION 

During inflation the vacuum fluctuation of each light (practically massless) field is promoted to a classical pertur- 
bation, on those scales which leave the horizon. Wc shall be interested in single-component inflation, during which 



1 The original PBH abundance calculations and early numerical studies, e.g. used the density contrast (on comoving velocity 

orthogonal slicings) to study the gravitational collapse. In a recent study [g (see also 0) a metric perturbation, the curvature 
perturbation on uniform density hypersurfaces was used instead, giving excellent agreement with detailed numerical investigations 0- 



the only relevant light field is the slowly-rolling inflaton field <fi. The perturbation of the field, 8<j), is related to the 
curvature perturbation 1Z on comoving, or uniform field hypersurfaces, by 

K = -§5<t>, (1) 
9 

where the right hand side of Eq. has to be evaluated on flat slices. In the leading order of the slow-roll approximation 
the field equation for S(j) on flat slices lives in unperturbed spacetime. (In other words, we can ignore the back-reaction 
of the metric perturbation.) 

In this section we recall the standard description of the behaviour during inflation of the vacuum fluctuation of </>, 
which actually applies to the vacuum fluctuation of any practically massless field. We use conformal time t to work 
with tp = aS(f>, where a is the scale factor. We take inflation to be exponential corresponding to constant Hubble 
parameter H, and take 

T = -(aH)-\ (2) 

so that r = corresponds to the infinite future. 

We consider the Fourier component of tp with comoving wave vector k, whose classical field equation is 

*£2i + (*-£)wM-o. (3) 

Well before horizon exit Eq. J3J reduces to 

^ + ^(k,r)=0. (4) 

This is the equation of a unit-mass harmonic oscillator with angular frequency k. 

Now we quantise the perturbation Stp while keeping the unperturbed field classical. At the quantum level, each 
observable corresponds to a hermitian operator. We use a hat to denote operators, so that the operator corresponding 
to <5y(x, t) is <5</?(x, r). The 'reality condition' ensuring that Stp(x, t) is hermitian is 

0(-k,T) = 0t(k, r ), (5) 

where the dagger denotes the hermitian conjugate. 

As is usual in quantum field theory, we adopt the Heisenberg picture, in which the state vector describing the 
system is time-independent and the operators corresponding to observables have the classical time-dependence. Then 
tp(k, t) is a solution of Eq. ®, and to satisfy the reality condition it must be of the form 

0(k, r) = tp(k, r)a(k) + tp*(k, r)at(-k), (6) 

where tp(k,r) is a solution of Eq. J2J), and a is a time-independent operator. 

To define tp(k, r) and a, we focus on the era well before horizon exit, and consider a region of space-time much 
smaller than H~ x (i.e. with At and Ax both much less than (aH)^ 1 ). In such a region space-time is practically flat 
(in other words, gravity is negligible) and we can assume that flat space-time field theory applies. The scale factor 
a is practically constant, and we can make the choice a = 1 instead of the usual do = 1. Then k is the physical 
wave- vector, and tp is the original field <p. 

Flat space-time field theory is constructed in such a way as to arrive at the particle concept. It boils down for scalar 
fields to noticing that each Fourier component has the dynamics of a unit-mass harmonic oscillator with position 
q = tp(k). The vacuum corresponds to all of the oscillators being in the ground state, and the equally-spaced excited 
energy levels of each oscillator correspond to the presence of 1, 2, • • • massless particles, each with momentum k. All 
of this is encoded by the flat space-time field theory mode function 



the commutation relations 

[a(k),at(k')]=^ 3 (k-k') [a(k),a(k')] = 0, (8) 
and the assertion that a annihilates the vacuum state, 



a (k) | vacuum) =0. (9) 



The solution of Eq. satisfying the initial condition Q is 



p(fc,r) = -Lf (l + ifcr)^e- ifcT . (10) 

We are supposing that the time-independent state vector of the Universe corresponds to the vacuum. At any 
instant, the vacuum is a superposition of states with definite values for the Fourier components y(k). Consider now 
some particular instant, and imagine that the Fourier components are measured. As the Fourier components are 
uncorrelated, the perturbation is gaussian so that all of its stochastic properties are determined by the two-point 
correlator, 

^(k)^(k')) = |^(fc)| 2 <5 3 (k + k'). (11) 

The spectrum of <p is defined as 

^ = 7T^ 2 ( fc )l> (12) 



2vr 2 a 2 



the normalisation chosen so that the mean square is 



m 2 ) = j Q n(k) T . as) 

Well after horizon exit it reduces to the famous result V$ = (H/2tt) 2 . 

The previous paragraph holds equally before and after horizon exit. The issue of quantum versus classical arises 
when we consider the evolution after the measurement has been made. The perturbation ip(\s) can be regarded as 
classical if there exist states, analogous to the "wave packets" that one considers in ordinary quantum mechanics, 
which are practically eigenstates of <£(k) over a long period of time. As we are working in the Heisenbcrg picture, this 
means that the perturbation is classical if and only if its operator has almost trivial time-dependence, 

0(k,r)~f(T)0(k,T Q ), (14) 

where / is a number and To is some initial time. Such is the case well after horizon exit, when |l4lll5| 

^(k,T)c^=L-[a(k)-at(-k)] , (15) 

but not earlier. 

This is the usual justification for regarding the field perturbation as classical well after horizon exit. It is not 
completely satisfactory since in the context of cosmology one can hardly talk about a complete measurement, but it 
works in practise. (Some of the conceptual issues are addressed in the extensive literature on decoherence, e.g. |lfij |. 
and on the quantum-to-classical transition |l7lllSj .) Our concern here is with those Fourier components <^(k), which 
at the end of inflation are still inside the horizon or at least not very far outside it. 

Limiting ourselves to the issue of black hole formation, we propose that the formation of the black holes should 
itself be regarded as a measurement. Accordingly, we will go ahead and calculate the mass function of the black holes, 
using the spectrum (|12|) exactly as if it referred to a classical quantity. In the end we will confine ourselves to the 
sub- horizon regime k 3> k en ^ in order to arrive at a simple and fairly model- independent result, but it is clear that 
within the context of a given model our methods would allow a continuous mass function to be calculated covering 
also the regime k ~ /c on d- 

Two notes of caution on the applicability of our results are in order here. Firstly, during the inflationary stage, 
since the amplitude of 1Z increases as k/aH increases, there is a scale below which it becomes greater than unity. It is 
given by k/aH = 1/IZhc, where IZhc is the value of 1Z at the horizon scale. Since we do not know how to deal with 
this scale, we confine our arguments to scales greater than this critical scale. Secondly, inside the horizon, the metric 
perturbations due to second order matter perturbations can in general become important. In the case of the quantum 
fluctuations, this seems to give the condition that the scale should be greater than V'((f>)/H 2 , or k/aH < V'/H 3 , 
where V is the inflaton potential and a dash denotes differentiation with respect to <j). Below this scale, the second 
order perturbations dominate. Hence the arguments in this article apply only to scales larger than these lower bounds. 



III. COSMOLOGICAL PERTURBATIONS 



A. The evolution of the Bardeen potential and the curvature perturbation 



The line element of the unperturbed spatially flat FRW universe is defined as 

ds 2 = a 2 (r)(-dr 2 + S lj dx l dx J ) , (16) 

where a(r) is the scale factor, r is the conformal time, Sij is the metric of the flat background, and x are the spatial 
coordinates. The perturbed metric in the longitudinal gauge, considering only scalar perturbations, is given by 

ds 2 = a 2 (t) [- (1 + 2*) dr 2 + (1 + 2$) 8 i:j dx l dx j } . (17) 

In this gauge the metric perturbations coincide with the gauge invariant Bardeen potentials defined in [l9| , W = $a 
and $ = &h ■ is the lapse function that determines the proper time intervals between zero shear hypersurfaces (the 
slicing of the longitudinal gauge) . $ is the curvature perturbation on slices of zero shear [23] • 

The relation of these two potentials to the density perturbation is very insightful. In the cases of a scalar field and a 
perfect fluid there is no anisotropic stress, and then <!> = — 'J. The potential if? is then related the density perturbation 
in the comoving gauge [l| 2 



p 3 \aH 

where the perturbation 8 k is the density contrast defined in the comoving slicing, p is the density of the unperturbed 
background, and H is the corresponding Hubble parameter, H 2 = 8nGp/3. Expression (|18fl is analogous to the 
Poisson equation in Newtonian gravity. 

The Fourier components of 1Z, defined in Eq. (JJJ above, are simply related to the Fourier components of the field 
fluctuations on flat slices 6<fr by 

n k = -^5cf> k . (19) 

The curvature perturbation lZ k sets to a constant value a few Hubble times after horizon crossing, getting a very 
small scale dependence determined by the value that H and <fi have at that time. Note that this result is obtained by 
considering the evolution equation for lZ k , but not from the approximate solution for ip k — ad<pk given by Eq. IjlOJI . 
On scales inside the horizon, the time dependence of lZ k is mainly given by the vacuum fluctuations in Ijl0|l . Then 
one may write 



1 h a 



K k (t) = { _ (20) 



/2fc3 



(fL iork«aH, 



where t* indicates the time when a given scale k leaves the horizon. 

The other quantity we are interested to determine from the vacuum fluctuations is the Bardeen potential ^ . In the 
absence of anisotropic stress the evolution of \& is given by the following first order equation 

2 -H~^ k + = -(l + w)K k . (21) 

where w = p/p. On super-horizon scales !Z k is constant, and the decaying mode can be neglected. Assuming w is 
only weekly time-dependent, the solution is given by 

*, = -^W (22) 
5 + 6w 



2 Here and in the following, quantities with specified wavenumber k are mode functions governing the evolution of the corresponding 
operators. Thus &k f° r instance is defined through IH1 with <p(k, r) replaced by 5(k, r) and <p(Je,r) replaced by S(k,r) = 5 k . 



During slow-roll inflation w ~ — 1 and therefore \^k\ <C \TZk\- Inside the horizon TZk varies with time, and the time 
derivative under slow-roll conditions from (|20|l is given by 

• 1H 3 ( k \ 2 f.k\ 

K k = ^=T 777 exp .— . (23) 



On the other hand, TZk can be written as pj 

where SPk is the pressure perturbation in comoving hypersurfaces, and SPk = Spk for the inflaton field. Comparing 
both expressions for TZk, the mode function for the Bardeen equation yields 

■MH 2 {l + w) ( . k \ 
= TT L^; ex P hrr? > (25) 



which together with the previous result on super-horizon scales, Eq. I|22[l . shows that ^>k during slow-roll is practically 
zero on all scales, provided that H 2 /\4>\ <C 1 which is required by observational consistency. 



B. Junction Conditions 



In this section we are setting the initial conditions for the radiation dominated epoch after the end of inflation. If 
reheating occurs almost instantaneously, then we can use the values of for TZk and VPfc at the end of inflation as the 
initial conditions for the radiation epoch, as we will show below. A hybrid inflationary model, in which the decay 
and thermalisation of the particles produced during reheating takes place very rapidly, is most appropriate to realise 
the definite values of these quantities. In these models, the end of inflation is determined by the critical value of the 
inflaton field <f) c , which triggers the rapid fall of the secondary field a to its true vacuum. Then the matching of the 
solutions must be done on comoving hypersurfaces, (^ c (x, i) = constant. 

The junction conditions require the intrinsic metric and the extrinsic curvature tensor of the comoving 3- 
hypersurfaces on which we match to be continuous [2lL I22I l23l | . From the requirement on the intrinsic metric we 
therefore find that TZ has to be continuous. The trace of the extrinsic curvature of a comoving hypersurface at a 
given point x, is the locally defined Hubble parameter if(x, t) = H(t) + 5H(x.,t). Therefore the perturbation SHk is 
continuous, and it can be written as 

87rG„ 2 /fc x 2 



2H6H k = — 5 Pk - - - K k , (26) 
3 3 \a J 

showing that the Bardeen potential ^fc, related to the comoving energy density 8pk by (|18|l . is also continuous. 

Since TZk and are continuous at the transition time in the comoving gauge, from Eq. H21[) ^fc must be discon- 
tinuous on that hypersurface. Its value can be determined at the beginning of the radiation epoch from J2U, setting 
^fe = and taking TZk given in 12Ufl . and is given by 

f + w )H (te) TZk (* e ) for k > a c H c , 

(te) Irad = < (27) 

[ -|(1 + w )H (t c ) TZk (**) for k < a c H c , 

where wq — p/p = c 2 s = 1/3 is the velocity of sound squared during radiation domination, a c and H e are the scale 
factor and the Hubble parameter at the end of inflation, respectively, and TZk (t* ) denotes the value of the curvature 
perturbation when the scale leaves the horizon given in (|20|l . 

To smoothly match the background scale factor between the inflationary stage and the radiation-dominated stage, 
it is necessary to make a constant time shift of the conformal time r. We set 

Hi^j for t <C t c , 
<>(t) = { ^ (28) 
lk forr>r c , 



where r e is the conformal time at the end of inflation. Thus, t in all the previous formulae should be replaced by 
r - 2r c . 



C. The Bardeen potential during radiation domination 



For a radiation dominated epoch the Bardeen potential is 



V k (r) 



ai (x cos x — sin x) + 02 (x sin x + cos a;) 



where a\ and a2 are integrating constants, and x = c s kr. 

From = one gets the following relation between a% and 02 



ai 



where x c = c s kr c , and A (fc, r e ) is given by 



A(k,T e ) 



-A(k,T c )a 2 , 

(x e sinx c + cosx e ) 



(x e cos x c — sin x Q 
Differentiating and evaluating Eq. (|29|) at r c = l/a c H c yields 



dr 



«2 



c s kTl 



cosx c + A (fc, t c ) sin x e 



Then taking the derivative of Eq. (|27[) . the integrating constants are 

cos x c + A (fc, r G ) sin x e 
(1 + w) 7?.fc (x c cos x c — sin x c ) , 



a 2 = c s krl 
3 



-- (1 + w) 

2 T n 



«1 



(l + ui)7?.fe (x e sin x c + cos x ) . 



Substituting the value of a\ and a 2 into Eq. (|29[) . one finds that 'J'fc during radiation domination is 

3(i + w )n k 



*fe (r) 



2a; 3 



(x — x c ) cos(x — x c ) — (1 + xx c ) sin(x — x c ) 



(29) 



(30) 



(31) 



(32) 



(33) 
(34) 

(35) 



The value of TZk in this last expression is given by Eq. (|20J) evaluated at r c . On scales k <C a c H c , vWlZk ~ (H 2 /<p) T=Tt , 
and it has a slight scale dependence, while for k > a e H e it has the same fixed value for all scales inside the horizon 
at r c . 

Long after the end of inflation, r 3> t c , and on scales bigger than the sound horizon, x — kc s r < 1, we recover the 
usual long wavelength result for given by Eq. (|22|l . Using wo = 1/3 and expanding the sine and cosine functions, 
we get from Eq. (|3"5jl 



(t) 



X'' 



y + X 2 X - O (X 5 ) 



l + ^ + 0(x 2 ) . 

which agrees to first order in x with the value obtained from Eq. 1221) setting w = 1/3. 



(36) 



IV. BLACK HOLE FORMATION INSIDE THE HORIZON AT THE END OF INFLATION 



A. Motivation 



The solution of the Bardeen potential for radiation domination, Eq. Ij35(l , shows that the oscillations of the potential 
are damped with time, and the maximum value that ^^(t) can reach in the linear regime is during the first oscillation. 
Previous studies indicate that the criterion for the BH formation from super-horizon scale perturbations is ^^(t) — 



0(1) at horizon re-entry. In fact, we expect this criterion to hold also for sub-horizon scale perturbations. The reason 
is as follows. Consider a perturbation localised on a physical length scale L. Eq. H18|) indicates 



y^GSptf-^i, (37) 

where SM is the total mass excess of the perturbation. On scales sufficiently smaller than the horizon the Hubble 
expansion can be neglected, and the condition 2 G SM > L is equivalent to the usual condition for the BH formation 
in the asymptotically flat spacetime. Therefore on all scales inside the horizon one expects that the condition \P ~ 1/2 
ensures that the perturbations have entered the relativistically nonlinear regime and started collapsing. On the other 
hand from Eq. (|18|l . ^f c = 0.5 corresponds to a value of the critical density contrast at horizon crossing of 5 C — 1/3, 
and this value of S c is slightly favoured from superhorizon PBH collapse numerical estimations 0, 0, 01 ■ 

One might wonder, on scales well inside the horizon if the usual approach, considering Jeans stability criterion for 
the perturbations in the linear regime, precluded black hole formation. However using the same criterion as above for 
horizon scale perturbations reveals that for W ~ 1/2 that would not be the case. 

B. Mass Spectrum Estimation 

The derivation of the mass spectrum of PBHs on scales smaller than the horizon at the end of inflation in Press- 
Schechter theory resembles the standard one where perturbations reenter the horizon after the end of inflation 0, [|| . 
Here we use the threshold criterion for PBH formation \& c = 0.5, and therefore to apply the Press-Schechter formalism 
we use the Bardeen Potential instead of the density field. 

The fraction of space of PBHs with mass larger than to for a Gaussian probability distribution is given by 



,oo : 

F (m, ^f c ) — / exp . 

Jq, c V27TCT* (to) \ 2cr^ (to 



* 2 



, (38) 



where a^{ni) is the mean square deviation of the Bardeen potential after smoothing the field on a given mass to. 
Using a top hat window function in Fourier space, the mean square deviation is approximately given by 

o\ (to) ~ 7>* (to) . (39) 

On sub-horizon scales the oscillating power spectrum Eq. 1)35(1 is damped out with time. However, the \&k ~ 1 
condition corresponds to strong gravitational effects and therefore to nonlinear evolution, and as a consequence, once 
a perturbation on a given scale reaches that value it evolves nonlinearly. For that reason, to estimate the number 
of PBHs formed on a certain scale at the end of inflation using Press-Schechter theory, one needs the mean square 
deviation of the smoothed field at approximately the time when the oscillating potential \I/k reaches the first maximum. 
The time dependence in this calculation is handled in a similar way as in the standard scenario. In that case as a 
perturbation enters the horizon black holes are formed if the amplitude of the perturbation is larger than the threshold 
value, and the size of the BH is determined by the time (and the horizon size at that time) at which it is formed. 
Here the size of a PBH is determined by the physical wavelength at which iff reaches the threshold value. 

In the short wavelength limit, x c 3> 1, the Bardeen potential in Eq. i|35f) can be written as 

2#e ( xl \ . . fix, 
. x sin (x — x c ) exp — 

In this limit, the Bardeen potential reaches the first maximum at a time r*, such that 



vf k (T)~^=%.(-^)sin( a; -x c )exp(^) . (40) 



7T 

(X - X c ) = fcc s (T* - T c ) = -. (41) 

As k approaches the value of horizon at the end of inflation k = a e H e , the maximum of the Bardeen potential cannot 
be calculated analytically since it involves a transcendental equation. Nevertheless, it can be shown that in that limit, 

(a;.-X e ) = 0(1). 3 



3 The difference (x* — x e ) approaches very rapidly to tt/2 as k increases. On scales between the causal and the sound horizon at the end 
of inflation, (x„ — x c ) < 1, while for smaller scales 1 < (x* — x c ) < w/2. 



The time at which the Vf^ potential reaches the first maximum determines the mass of the PBH, m, and the 
smoothing scale for the field. The mass of the PBH is given then by 



47r /a*\ 3 4-7T fa c 



Xo. + 1 



(42) 



where is the energy density of the unperturbed background at r = r*, and (a^/ac) 1 / 2 ~ 1 + l/x e . 

To estimate the mean square deviation at the maximum, first consider the time dependent expression smoothed on 
the mass scale to, which from Eqs. I|35|) and (|39[) is given by 



a% (m, t) 



iVn (x c ) 



(x — x c ) cos(x — x c ) — (1 + xx c ) sm(x — x c ) 



(43) 



where Vfc(x e ) is the power spectrum of the curvature perturbation, defined by V-n{k,T) — k 3 /2n 2 < |7?.fc(r)| 2 >, 
evaluated at the end of inflation. From Eq. Q20JI. V-n(xc) can be written for x c > 1 as 



Xr 

1 + — 



(44) 



Evaluating the mean square at the approximate maximum, (x* — a; c ) ~ 1, Eq. 143|l for x c > 1 is approximately given 
by 4 



ay(m) ~ 2V^ 2 {x c ) 



2A 1 ' 2 I 1 + 



7? 



2\ 1/2 



where *4^ 2 = {H 2 / <j) c ) /2ir. Note that the value of converges in the limit x e 3> 1 to give 



ay = 2^ /2 /c s = 2v / 3.A^ 



1/2 



(45) 
(46) 

(47) 



This is a very important feature for the formation of PBHs inside the horizon, as the number of mean deviations 
Ny = ^ c /ay reaches a minimum value. 

Once we have determined ty c and ay , we can write the mass fraction of PBHs per logarithmic mass interval 



dF (to, # c 
dm 



2 mday 
7r ex 2 , dm 



■ exp 



* 2 

_ C 

'24 



(48) 



The derivative of the mass with respect to the comoving scale is 

dm 47r /a c \ 3 1 (2 + 3:r e ) 



dm 47T /a c \ J 1 (^z + o^cj 
~dk ~ 3~ \k) k x f 1 + JL 



to 2 + 3x e 



fc 1 + x 

On the other hand the derivative of the mean square deviation with respect to the comoving wavenumber yields 

day day dx c ay 6x 2 — x + 1 
dk dx e dk k (1 + x 2 /c 2 ) (1 + x e ) 

These two last expressions allow us to write the derivative of the mean square deviation in the following form 

day day dx e dk ay 



(49) 



(50) 



dm dx e dk dm 



F{x e ) , 



(51) 



4 This approximation breaks down for x c < 1, but the mean square value converges very fast to the value given in Eq. 1451 as x c 
increases over the value 1. 



where T{x e ) is defined as 



6x 2 — x c + 1 
(1 + x 2 /c 2 ) (2 + 3z e ) 



, a,; , — ■ (52) 



Taking expression 1)51(1 the mass spectrum can be written as 

dn (to, go) [2 p ( N 2 \ 

where iV* = ^c/cr* is the number of mean deviations of the Gaussian distribution of the metric potential and p 
is the energy density of the background at formation time. 

V. DISCUSSION AND CONCLUSION 

The discussion of the results becomes more transparent if we introduce the quantity (m 2 / p)dn/dm, the mass fraction 
of PBHs per logarithmic interval of mass. Unlike the number density dn/dm, the mass fraction is finite on all mass 
scales and therefore the total mass of PBHs is finite as expected. From Eq. 1(53(1 the mass fraction (m 2 / p)dn/dm is 

m 2 dn(m) . . 1T / N 2 \ 

= \\-T (x e ) JV* exp -- ^ , (54) 



p dm V 7r V 2 

valid for A-R 3 ^ 2 m c < to < 0.1to c , where m c is the horizon mass at the end of inflation, m = (47r/3)/9c-ff<r 3 - The lower 
bound is the required limit according to cosmological linear theory. The upper bound on the validity of Eq. 1(54(1 . 
to < 0.1to c , stems from the following requirements. The relation between x e and the mass of the black hole inside 
the horizon, m, is given by Eq. ((42(1 . and we get 

x Ux c + l)^c^. (55) 
to 

Then the mass fraction inside the horizon in Eq. 1(54(1 is valid for x e > 1 or to < 0.1m o . For simplicity, in the 
following plots we have taken x e = c s (m/m ) -1 / 3 . This amounts to neglecting the term ln(l + l/x c ) when we take 
the logarithm of m/m e , and again the approximation breaks down for x e < 1. 

Using the upper bound on the mean deviation of the density contrast at horizon crossing in 0,0] , &hc < 0.04, 
we can set an upper bound on the amplitude of the fluctuations on sub-horizon scales A-r 1 ^ 2 — (9/4)cthc J$ 0.09. The 
following graphs are shown highlighting the regime in which our calculation is valid. The mass fraction (to 2 / p)dnj dm 
inside the horizon increases rapidly with the exponential factor as the number of mean deviations N\j, decrease. Hence 
the mass fraction rises sharply with k, as shown in Fig. ^ with a solid line for a threshold value of $> c = 0.5 and 
an amplitude A-r 1 ^ 2 = 0.09. In Figs. El and 01 we show the logarithm of the mass fraction for a threshold value of 
ty c = 0.5 and an amplitude A-r 1 ^ 2 = 0.09, and also a more conservative amplitudes of A-r 1 ^ 2 = 0.05 and 0.02. The 
extension of the mass fraction in equation 1(54(1 beyond the limit to = A-R 3 ^ 2 m c is shown in dotted lines in Figs.^and 
OH reaching a maximum and eventually converging to zero. 

There are several things worth noting in Figs. |21 and |3 (i) The maximum is displaced to smaller scales as the 
number of mean deviations increases, (ii) As increases the mass fraction maximum value falls dramatically. 
The curve for Ar = 0.01, which is not shown in Fig.|2 has a maximum with an approximate value of 10 -60 . (hi) The 
maximum of the mass fraction and the lower bound nearly coincide for A-r 1 ^ 2 — 0.09 on a mass scale of approximately 
to ~ 7 x 10~ 4 to o . However, Fig. [21 shows that the drift of the maximum into smaller scales is more rapid than the one 
of the current lower bound as A-r 1 ^ 2 decreases its value. 

An estimate of the total mass fraction inside the horizon is provided by integrating Eq. ((54(1 over the relevant range 
of masses. Considering a sharp cut-off at m = A-R 3 ^ 2 m c , the total mass fraction gives an upper limit of 4 x 10~ 2 for 
a threshold value ^ c = 0.5 with an amplitude of A-r = 0.09. If we extend the validity of Eq. ((541) into the ultraviolet 
regime, the upper limit of the total mass fraction increases to approximately 0.14. The difference between these two 
estimates of the total mass fraction increases considerably as A-r 1 ^ 2 decreases. Therefore, the behaviour of the mass 
function in the ultraviolet regime needs further investigation. We postpone this issue and the detailed calculation of 
the observational effects of sub-horizon PBH formation, in particular the implications for PBH abundance calculations, 
to a future publication. 
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FIG. 1: The mass fraction (m 2 / p)dn/ dm versus y = ln(m/m c ) for ^ c = 0.5 and A-r 1 ^ 2 = 0.09. The dotted line correspond to 
the extension of the mass fraction in Eq. 154H beyond the lower bound y ~ 7. 




FIG. 2: \og 10 ({m 2 /p)dn/dm) versus y = ln(m/m e ) for $ c = 0.5. The dotted line corresponds to the amplitude Atz 1 ^ 2 = 0.02, 
the dashed line to Aiz 1/2 = 0.05, and the full line to Aiz 1/2 = 0.09. 




FIG. 3: log 10 ((m / p)dn/dm) versus y = ln(m/m e ) for * c = 0.5 and Ar 1 = 0.02, 0.05 and 0.09. Here the dotted lines 
correspond to the extension of the mass fraction in Eq. Q54| beyond the lower bound m = Aiz'^m^. 
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